Quantum escape of the phase in a strongly driven Josephson junction 
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A quantum mechanical analysis of the Josephson phase escape in the presence of both dc and ac 
bias currents is presented. We find that the potential barrier for the escape of the phase is effectively 
suppressed as the resonant condition occurs, i.e. when the frequency uj of the ac bias matches the 
Josephson junction energy level separation. This effect manifests itself by a pronounced drop in the 
dependence of the switching current I s on the power W of the applied microwave radiation and by 
a peculiar double-peak structure in the switching current distribution P(I a ). The developed theory 
is in a good accord with an experiment which we also report in this paper. The obtained features 
can be used to characterize certain aspects of the quantum-mechanical behavior of the Josephson 
phase, such as the energy level quantization, the Rabi frequency of coherent oscillations and the 
effect of damping. 
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Great attention has been devoted to the experimental 
and theoretical study of Josephson junctions in the pres- 
ence of externally applied microwave radiation. The well 
known effects are Shapiro steps in the current-voltage 
characteristics, and features due to photon-assisted tun- 
neling of quasiparticles Moreover, microwave radi- 
ation can be used to probe the macroscopic quantum- 
mechanical behavior of the Josephson phase such as the 
energy level quantization 0, m and the recently ob- 
served quantum coherent Rabi oscillations [a, Kg. The 
observation of these effects has proven that the laws of 
quantum mechanics can be applied to dissipative sys- 
tems. The interest in this field has also been boosted by 
the possibility to use Josephson junctions as supercon- 
ducting qubits 

The quantum mechanical properties of the Josephson 
phase can be studied most directly when the Josephson 
phase interacts resonantly with the external microwave 
radiation la, K* L3 • Such an interaction leads to the 
resonant absorption of photons, which induces transitions 
between the energy levels of the Josephson phase and 
thus to a pronounced increase of the population of excited 
levels. 

It is also well known that thermal and quantum fluctu- 
ations lead to a premature switching of the junction from 
the superconducting state to the resistive state at a ran- 
dom value of dc bias I s < I c , where I c is the critical cur- 
rent in the absence of fluctuations 0,12)111. I* 1 this case, 
the escape of the Josephson phase can be characterized 
by the switching current distribution P{I S ) measured by 
ramping up the current at a constant rate. The resonant 
absorption of photons in a Josephson junction results in 
an enhancement of the rate of the phase escape at a par- 
ticular dc bias I r 0j S S 0| • The value I r is found from 
the resonant condition as the frequency uj (or nuj, where 
n is an integer number, for multiphoton absorption Q) 
matches the energy level separation 



Here, E n ^ m (I r ) are the quantized energy levels of the 
Josephson phase in the presence of an externally applied 
dc bias I r . The enhancement of the escape rate manifests 
itself by additional microwave-induced peaks in P(7»). 
This effect has been analyzed in detail in Refs. 

The crucial condition to observe such an enhancement 
of the escape rate is the small-signal limit, in which the 
microwave radiation only changes the population of indi- 
vidual levels in the well. To analyze this condition quan- 
titatively we note that the fluctuation induced escape of 
the Josephson phase from excited levels has to be large. 
In the quantum regime the escape probability is governed 
by the exponent of the quantum tunneling rate which is 
proportional to the ratio of the barrier height to the at- 
tempt frequency 

£j(l-7) 3/2 

where Ej and uj p are the Josephson energy and the 
plasma frequency of the junction, and 7 = I s / 'I c . To 
observe fluctuation-induced escape of the phase this ex- 
ponent has to be of the order of unity, i.e. 1 — 7 < 
(hujp/Ej) 4 '' 5 . At the same time to observe the resonance 
the drive frequency uj has to fulfill the resonance condi- 
tion a; = uj p [2(1 — 7)] 1 ^ 4 . Thus, an enhancement of the 
escape that can be directly interpreted as a microwave- 
induced increase of the excited level population, occurs 
approximately as 



< 
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(2) 



KuJ = E n (I r ) ~ E m (I r ) 



(1) 



As we turn to a larger microwave frequency or to 
smaller value of the critical current density, the condi- 
tion J2j is not valid any more, and a simple increase of 
the population of excited levels can not explain the res- 
onantly induced escape of the Josephson phase. In this 
paper we show that in the limit (uj/uj p ) 5 3> (huj p )/Ej, 



2 



the presence of microwaves leads to an effective suppres- 
sion of the potential barrier for the escape of the phase. 
This suppression is particularly strong when the reso- 
nance condition is satisfied. We show that this resonant 
effect does not only lead to a double peak structure of the 
P(I 8 ) distribution at the resonance but is also accompa- 
nied by a sharp drop in the dependence of the average 
switching current (I S (W)) on the power W of the applied 
microwave radiation. Note here, that the analysis of this 
effect based on the classical nonlinear dynamics for the 
Josephson phase has been presented in [8j . 

To analyze the problem quantitatively we write down 
the time-dependent potential energy of the Josephson 
phase 4>(t) in a small Josephson junction in the presence 
of both dc bias current / and the ac bias with amplitude 
7/ induced by microwaves: 



U (<f>)-Ej^-wx{Ljt)<f> 



Ej ( cos0 + 



(3) 



Here, Uo((f)) is the potential energy of the Josephson junc- 
tion in the absence of ac bias. 

In analogy to the classical nonlinear problem |g, we 
write the coordinate 4>{f) as a sum of two terms: a quickly 
oscillating resonant term and a term <po(t) slowly 
varying in time. Next, the quantum-mechanical average 
of (£(£)) is calculated using perturbation theory with re- 
spect to the small amplitude 77. Moreover, we truncate 
the dynamics of the Josephson phase in the potential 
Uo(<p) to the two energy levels E n ^ m which resonantly 
interact with the ac bias. We take into account the 
damping parameter a that determines the strength of the 
resonant interaction. This damping can be governed by 
the quasi-particle resistance of the junction itself and/or 
dominated by the impedance of the bias leads at the fre- 
quency io [3IIH. lil Ho| . In linear approximation we obtain: 
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, (4) 



where the energy level differences E nm (I) = E n (I) — 
E m (I) and the matrix elements f nm = (n\<j)\m) depend 
on the dc bias current. The values of E n (I) and the 
corresponding eigenfunctions \n) are found as solutions 
of the Schrodinger equation 



2Ej 



(5) 



Substituting the term in the expression for the 

potential energy U ((f)) and taking the average over time 



we obtain the effective potential energy U e ff{<j>o) m the 
form 



U eff (M = ~Ej 



— + COS 0o (1- 
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The dependence of the switching current I S (W) on the 
microwave power W = krj 2 /2, where k is the microwave 
coupling coefficient is obtained by making use of the con- 
dition that E/eff (</•()) has no extreme points. Thus, the 
shift in switching current SI S (W) — (I c — I S (W))/I C is 
determined by a solution of the transcendental equation 
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Note here, that as we derive Eq. (JJJ all fluctuation effects 
are neglected. 

A good approximation of the dependence i5/ s ( W) can 
be obtained in the harmonic oscillator model. In this 
case, the term with (n = 0, m = 1) in right-hand part 
of the Eq. (JJJ is most important, and the energy level 
difference is Eqi = hujp(2SIs) 1 ^ 4 ■ The typical calculated 
dependence of the switching current on the microwave 
power is presented in Fig. ^ The most peculiar feature 
of the dependence I S (W) is a sharp drop as the microwave 
power is close to the critical value W cr . The appearance 
of the critical value of microwave power has a simple 
physical meaning. As the microwave power W is less 
than W cr the resonant interaction is a weak, and the ac 
induced suppression of the barrier is not enough to allow 
for the Josephson phase to escape. The I S {W) curves also 
display a weak dependence on W in the limit W < W cr , 
and this dependence becomes stronger in the limit of rel- 
atively high microwave power W > W cr - Moreover, due 
to the multi- valued character of I s (W) there is a particu- 
lar region of microwave power where a switching current 
distribution P(I S ) with two peaks should be observed. 
Notice that, while the width of the peak corresponding 
to the larger value of the switching current / s is due to 
the presence of fluctuations (thermal or quantum), the 
width of the lower peak is determined by the damping 
parameter a. The magnitude of the critical current drop 
becomes smaller as the damping parameter a increases, 
see Fig. Q 

The above analysis can be extended to include fluctu- 
ations. The presence of thermal fluctuations leads to a 
premature switching of the junction from the supercon- 
ducting state with respect to the value I S {W) obtained 
from Eq. J7J). Taking into account thermal fluctuations 
we derive a transcendent equation that is similar to (Q) 

(SI S (W)) = (SI s (0))+ 
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FIG. 1: Microwave power dependence of the switching cur- 
rent I S (W) of a small Josephson junction. Dashed lines for 
a — 0.06 bound the region in which switching current distri- 
butions with two peaks should be observed. The frequency 
u> of the ac drive is — 0.4 and the damping parameter is 
varied from a — 0.04 lo p to a = 0.08 ui p . These correspond 
to the junction's effective quality factor Q = q~ 2 ui 2 ranging 
from 100 to 25. 
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where (5I S (0)) is the fluctuation-induced shift of the 
switching current in the absence of ac bias. For tem- 
peratures above the crossover to the quantum regime 
(SI s (0)) = [^ln(^f)} 2/3 , where / is the bias cur- 
rent ramp rate p}. By solving the Eq. J8J) for a — const 
we find that thermal fluctuations lead to smearing of the 
drop in the mean switching current and also shift the 
critical microwave power W cr to smaller values, as illus- 
trated in Fig. 

Equations {7} and JSJ can also be used in the more 
complex case of a strongly non harmonic potential well 
Uo((/>). In such a regime, other terms in Eq. J7J) become 
important. For example, by taking into account E m tran- 
sitions we should expect, according to Eq. Q, an addi- 
tional jum p (of much smaller value) in the dependence 
I S (W) 11]. The microwave power W required to observe 
this jump is less than W cr , the critical microwave power 
causing the main drop in the I S {W) curve. The multi- 
photon absorption of the Josephson phase observed in 
Ref. 0, can be also be treated in the framework of this 
analysis. In this case Eqs. J7J and JSJl can still be 
used with the substitution of to by mo. Instead of the lin- 
ear dependence of on the amplitude 77, one should 
expect the dependence (£(t)) ~ Jn{^^), where J„ is the 
Bessel function of order n |l2j ■ 

The coupling coefficient k can be found by comparison 
of the observed (I S (W)) dependence with the one calcu- 
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FIG. 2: Temperature dependence of the switching current of 
a small Josephson junction in the presence of ac bias at the 
frequency u = 0.4 u p . The reduction of the critical current 
due to thermal fluctuations is taken into account. The damp- 
ing parameter is fixed a = 0.04 u p . Temperatures shown are 
0.1 K, 0.2 K and 0.3 K. 



lated in Eq. iJSJ . It allows to calculate the crucial param- 
eter of coherent quantum-mechanical behavior, i.e. the 
microwave power dependent Rabi frequency of coherent 
oscillations as 



— (jJpfnm\ —j — 



2W 



(9) 



It is worth to note that our analysis is also applicable to 
the problem of the microwave-induced escape of a pinned 
magnetic fluxon in an annular long Josephson junction, 
the dynamics of which can be map ped to the dynamic of 
a small Josephson junction [illl3|. 

The analysis presented here was motivated by mea- 
surements of the switching current distribution P(J S ), 
which we performed on a small (5x5 fim 2 ) Nb/AlO^/Nb 
Josephson tunnel junction. This junction had a 
fluctuation-free critical current of I c = 278.45 /xA, an ef- 
fective capacitance of C = 1.61 pF and a quality factor of 
Q w 45 ± 5 J]. P(I S ) distributions have been measured 
at a bias current ramp rate of 0.245 A/s at the temper- 
ature T = 100 mK, which is well below the temperature 
T* sa 280 mK, corresponding to the crossover between 
thermal escape and quantum tunneling for this sample. 
Microwaves in the frequency range between 10 GHz and 
38 GHz have been applied to the sample and the mi- 
crowave power was varied over a wide range. In Ref.0 we 
have focused on the regime in which for small microwave 
powers the population of an excited state is enhanced and 
subsequent tunneling from the excited state is observed. 
In those experiments the drive frequency ui has been cho- 
sen to satisfy the condition ©. Here, n is the number 
of photons involved in the transition. In this case the 
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P(J S ) distributions were only weakly perturbed by the 
microwaves. 

Here we present experiments for relatively large mi- 
crowave frequencies nuj/uj p . To observe a microwave- 
induced resonance under this condition, a larger mi- 
crowave power has to be applied. In this case, the P{I S ) 
distributions are strongly perturbed by the microwaves. 
In Fig. the dependence of the most probable switching 
current I p on the microwave power W is plotted for a mi- 
crowave frequency of v = oj/(2tt) = 24 GHz. In Fig. [SJd, 
the P(I S ) distributions are plotted for the three selected 
microwave powers indicated in Fig. [3]}. It is observed that 
at small microwave powers the switching current distri- 
bution initially shifts to lower bias currents without de- 
veloping a resonant double peak structure. At larger mi- 
crowave powers - when the resonance condition is fulfilled 
- a double-peak distribution occurs. It is accompanied 
by a pronounced drop in critical current as the microwave 
power is further increased. The resonance presented here 
corresponds to n = 2 photons with 2lu/u> p — 0.41 [l4| . 

The observed features in the P(I S ) distribution can not 
be explained by a simple enhancement of the Joseph- 
son phase escape due to a microwave induced increase 
of the population of an excited state 0, 0, • We find 
a particular microwave power range in which two states 
(resonant and non resonant ones) coexist, indicating a 
microwave-induced bistability of the junction predicted 
by our theoretical analysis presented above. The charac- 
teristic behavior in the region of W ~ W cr is also in a 
good agreement with analysis above, compare Fig. and 
Fig. Hi 0. 

In conclusion we have carried out a quantum- 
mechanical analysis of the Josephson phase escape in 
a small Josephson junction in the presence of both dc 
bias and microwave currents. We find that in the ab- 
sence of fluctuations and for large microwave frequency 
the dependence I S (W) displays a pronounced drop at res- 
onance and thereafter decreases smoothly. This behavior 
is explained by the effective suppression of the poten- 
tial barrier by microwaves and is in a good accord with 
experiments carried out with small Josephson junctions. 
Similar behavior has been observed for a pinned Joseph- 
son vortex in a long annular Josephson junction 13] . 
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